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ABSTRACT 

We present new EVLA continuum observations at 7 mm of the 253-1536 binary disk system in the 
Orion Nebula Cluster. The measured fluxes were combined with data in the sub-mm to derive the 
millimeter spectral index of each individual disk component. We show how these observations can 
be used to test the models of dust evolution and early growth of solids in protoplanetary disks. Our 
analysis indicates that the disk with lower density and higher temperature hosts larger grains than 
the companion disk. This result is the opposite of what predicted by the dust evolution models. The 
models and observational results can be reconciled if the viscosity a-parameter differs by more than 
a factor of ten in the two disks, or if the distribution of solids in the disks is strongly affected by 
radial motions. This analysis can be applied to future high-angular resolution observations of young 
disks with EVLA and ALMA to provide even stronger observational constraints to the models of dust 
evolution in protoplanetary disks. 

Subject headings: circumstellar matter — protoplanetary disks — planets and satellites: formation — 
submillimeter: planetary systems 



1. INTRODUCTION 

Young circumstellar disks around pre-Main Sequence 
(PMS) stars are thought to host the birth of planets. 
According to the core accretion scenario, the formation 
of planets involves a variety of physical mechanisms from 
the coagulation of sub-/im sized particles up to the gas- 
accretion phase leading to the build up of gas giants. One 
of the most critical steps which is left to be explained is 
how this growth of solids proceeds from mm/cm-sized 
dust grains to km-sized rocks, called planetesimals. This 
is crucial because it is from the gravitational-driven col- 
lision of these planetesimals that both the rocky Earth- 
like planets and the rocky cores of giant planets are ul- 
timately supposed to be formed. Different mechanisms 
have been proposed in the literature to explain the for- 
mation of planetesimals. In general, these mechanisms 
induce the local accumulation of particles with sizes of 
~ 1 — 100 mm at densities which are high enough to 
make these clump s unstable to gravitational collapse 
(|Chiang fc Youdinll2010[ ). It is therefore very important 
to well characterize the early phases of coagulation of 
small dust grains in order to set the right initial con- 
ditions for the mechanisms governing the formation of 
planetesimals. 

Recently, sophisticated models of dust evolution in 
disks including several effects like dust coagulation, 
fragmentation and radial motions, have been built 



(jBrauer et al.l 120081 : iBirnstiel et al.ll2010at ). Continuum 
observations in the (sub-) millimeter constrain fundamen- 
tal parameters of the dust population in the disk, e.g. the 
total dust mass, the radial-dependent dust surface den- 
sity, and the size distrib ution of dust grains at sizes of 
about ~ 0.1 — 10 mm (see Willi ams fc C ieza 2011). These 
observational constraints can be used to test the models 
of dust evolution and s hed light onto the first stages of 
planetesimal formation (|Birnstiel et al.ll2010bl ). 

The recent initial upgrades of the Very Large Ar- 
ray into the Ex panded Very Large Array (EVLA; 
iPerlev et all 1201 11 1 allowed a significant increase in con- 
tinuum sensitivity at long wavelengths. We present 
new EVLA observations at about 7 mm of the 253- 
1536 binary system in the Orion Nebula Cluster. The 
projected angular separation of the two PMS stars is 
about 1.1 arcsec which corresponds to a projected phys- 
ical sepa ration of 460 A U at the ONC distance of ~ 
420 pc (jKraus et al. 2009). Optical images with the Hub- 
ble Space Telescope (HST) has revealed the presence of a 
large disk s een in absorbtion around the east companion, 
253-1536a. iMann fc Williamsl (|2009ft observed the binary 
at 0.88 mm with the Sub-Millimeter Array (SMA) and 
estimated for the 253-1536a disk a mass of about 0.07 M 
which makes this disk the most massive ever observed in 
the ONC. They also detected a fainter disk around the 
other PMS star, 253-1536b, and found a ratio for the two 
disk masses of about 4. 
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2. OBSERVATIONS 
2.1. EVLA 

We obtained new EVLA data under the project 1 OB- 102, 
with a total observing time of 6 hours. At the time of 
observations the EVLA was in the C configuration, and 
provided a total bandwidth of 256 MHz, that we centered 
at the frequency of 43.280 GHz (6.9 mm) in the Q-band. 

Passband and absolute flux calibration was performed 
through observations of the QSO 3C147 (J0542+4951). 
Amplitude and phase calibration was obtained by observ- 
ing the QSO J0607-0834. The measured 6.9 mm-flux 
density of J0607— 0834 range from about 1.9 Jy to 2.1 Jy 
during the three observing runs, thus showing a discrep- 
ancy in the fluxes which is within the ~ 10%-uncertainty 
on the absolute flux, as typically estimated for EVLA 
observations in the Q-band. 

The raw visibilities were calibrated using the Com- 
mon Astronomy Software Applications (CAS A) pack- 
age. Maps of the continuum emission were derived by 
adopting natural weighting to maximize the observations 
sensitivity, and photometry was obtained through Gaus- 
sian fitting in the image plane using the imf it task in 
CAS A. The resulting synthesized beam has sizes of about 
0.77" x 0.55". The measured rms- noise of the final map 
is about 35 /iJy. 

2.2. X-Shooter 

We observed 253-1536 with XShooter at the ESO-VLT 
on September 27, 2010. We obtained complete spectra 
from 330 through 2500 nm at a spectral resolution in 
the range 5000-9000. We aligned the - 11 x 1.0" arc- 
sec slit along the binary position angle to observe the 
two components simultaneously. We performed an on-slit 
dithering pattern with 8 exposures of 190s each, to avoid 
saturation on the brightest component and to efficiently 
remove sky emission, bad pixels and cosmic ray hits. We 
reduced the data using the standard XShooter pipeline 
recipes, and IRAF for extraction of the ID spectra and 
subsequent analysis. We used spectra of LTT 1020 and 
Hip 025028 (B3V) to correct for instrumental response 
and telluric absorption. The spectrum of the bright com- 
ponent is detected with very high signal-to-noise across 
the whole spectral range, while the faint component of 
the binary system is only detected at K and H bands and 
cannot be easily separated from the glow of the other 
component at shorter wavelengths. 

We detect numerous absorption and emission features 
in the two spectra. Nevertheless, for the purpose of this 
Letter we just analysed the spectra to attempt a spec- 
tral classification of the two components of the binary. 
We compa red our spectra with those from the spectral 
l ibrari es of I Allen fe Stroml (jl995f ) and IWallace fe Hinkld 
(1997). We find that the bright stellar component of 
the binary, i.e. 253-1536b, is consistent with previous 
classification as an M2 star both based on the optical 
and near infrared spectrum. The faint component, 253- 
1536a, does not show prominent absorption lines. In 
particular, it is laking CO, Fe, Al and Ca features that 
are commonly observed in K-band spectra of late type 
dwarfs. Overall the spectrum seems to be consistent 
with heavily veiled and extincted late F and G types. 
In this paper we will consider a spectral type of G2 for 
the 253-1536a PMS star. 



3. RESULTS 

The contours in Figure [2] show the EVLA map at 6.9 mm 
superimposed to the HST image in Ha. Both the disks 
are clearly detected at 6.9 mm, and the measured fluxes 
are listed in Table [TJ together with the ones measured 
with the SMA array at 0.88 mm ([Mann fe Williamsl 
2009). Continuum observations of young circumstellar 
disks in the sub-mm probe the emission from dust par- 
ticles. At longer millimeter wavelengths the dust emis- 
sion can be contaminated by ionized gas, with ioniza- 
tion driven by X-ray /UV-radiation generated either by 
near by young massive sta rs, e.g. in the Trapezium Clus- 
ter ([Williams et al.l 120051 ). or by the PMS star at the 
center of the disk, as seen for disks in low mass SFRs 
(jRodmann et al.ll20Q6D . 

In the case of the 253-1536 system, the distance to 
the 1 Ori C massive star, i.e. ~ 1 pc in projection (see 
Fig. [1]), is too large to make the gas in the disks being sig- 
nificantly ionized by the 1 Ori C UV-radiation. This is 
indicated by the lack, in these ONC peripheral regions, of 
photoevaporating proplyds with bright ionization cusps 
pointing; toward the dire ction of 6 1 Ori C in the HST 
images ([Ricci et al.l [2008). At the same time, a bright 
rim of ionized gas is visible in both the HST and EVLA 
maps on the east side of the binary system (Fig. [2]). The 
star which is causing this feature is probably NU Orionis, 
an early B-type star located toward the direction of the 
bright rim, which is al so responsable for the excitation 
of the M43 HII region ([Smith et al.ll2005D . 

To evaluate the relative impact of dust and ionized gas 
emission at a given frequency, the slope of the SED can 
be used: dust thermal emission presents a steep spec- 
trum, with values of the spectral index a mm (F v oc v a ) 
between about 2 and 3, whereas the ionized gas emission 
is characterized by a much more shallow dependence on 
frequency, e. g. a < 1 for a symme tric, ionized, optically 
thick wind ([Panagia fe Felli 1975), and a < for opti- 
cally thin emission ( Mez geT et al .1 Il967l ) . For the 253- 
1536 system, iRicci et all ([20111 ) measured a spectral in- 
dex of 2.6 ± 0.4 between 0.88 and 2.9 mm. Since the 
observations at 2.9 mm did not spatially resolve the bi- 
nary system, this value is descriptive of the whole system 
rather than of the two disks separately. By adding up the 
6.9 mm-flux densities for the two disks, the spectral index 
of the whole system is 2.0 ± 0.5 between 2.9 and 6.9 mm 
and 2.3 ± 0.2 between 0.88 and 6.9 mm. Both the mea- 
sured absolute values of the spectral indeces, and the fact 
that no statistically-significant changes are seen in the 
slope of the spectrum between 0.88 and 6.9 mm suggest 
that the dust dominates the total sub-mm/mm emission. 
This is consistent with what found also for 216-0939, an- 
other HST-identified d isk in the northern outskirts of the 
ONC ([Ricci et al.ll2Qllh . In the following we will there- 
fore assume that the measured sub-mm/mm fluxes are 
only due to thermal emission from dust particles in the 
disks. 

4. DISCUSSION 

4.1. Comparison of grain growth in the two disks 

Since the dust emission of disks at these long wavelengths 
is mostly optically thin, the spectral index a of the (sub- 
)mm SED is a proxy for the spectral index j3 of the dust 
opacity coefficient (n v oc i/P), which in turn carries the 
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TABLE 1 
Disk fluxes and spectral indeces. 



Source 


^6. 9mm 


rms 6 .9mm 


^0.88mm 


mis o.88mm 


Q!0.88-6.9mm 




(mJy) 


(mJy) 


(mJy) 


(mJy) 




253-1536a 


1.10 


0.035 


135.0 


1.0 


2.34±0.16 


253-1536b 


0.34 


0.035 


38.4 


1.0 


2.30±0.19 



Column 6: Spectral index between 0.88 and 6.9 mm; reported errors take into account a 10 %-uncertainty on the absolute flux scale at 
both 0.88 and 6.9 mm. 



information on grain growth in the disk. For example, 
in the Rayleigh- Jeans regime for a completely optically 
thin emission, f3 = a — 2. However, since the hypotheses 
of completely optically thin emission in the Rayleigh- 
Jeans regime are not always realized even at these long 
wavelengths, modelling of the disk emission is required. 

We modelled the measured disk sub-mm SED of 
the 253-1536 binary system by using a modified ver- 
sion of the two-layer pass i vely irradiated disk m odels 
(iChiang fe Goldreichl fl99^jDullemond et al.l l200lh . as 
described in iRicci et al The effective tem- 

perature of the two stellar companions were derived by 
converting the spectral t ypes adopted in Secti on [2^21 with 
the temperature sca l e of iLuhman et al.l ([20031 ) . By using 
the lPalla fc Stahlerl (|1999[ ) PMS evolutionary tracks and 
new multi-band op tical photometry and spectroscopy, 
iDa Rio et al.l (|2010l ) recently derived a new estimate for 
the ONC mean age of about 2 Myr. With these same 
PMS evolutionary tracks and age, the values for the stel- 
lar masses and luminosities are about 2 M , 12 L for 
253-1536a and about 0.3 M , 0.2 L for 253-1536b, re- 
spectively. 

The angular resolutions of the EVLA and SMA ob- 
servations do not allow us to properly constrain the ra- 
dial profile of the dust surface density in the two disks. 
For this reason, we consider in this analysis disks with 
truncated power-law surface densities (Edust oc r~ p for 
r < R out and Hdust = for r > R ou t) with possible 
p- values between and 1.5, as obtained through high- 
angula r resolution sub-mm imag ing of disks in nearby 
SFRs ([Andrews fc Williamsl 120071 ). As for the disk outer 
radius i? ut> the HST and SMA observations constrained 
a value of R out w 280 AU for the 253- 1536a disk and 
derived an upper limit of 60 AU for the unresolved 253- 
1536b disk. Since the R ou t of the 253-1536b disk is not 
determined by the observations, in the following anal- 
ysis we will consider two different possibilities, namely 
Rom =40 and 60 AU. Smaller disks with R out < 30 AU 
always fail to reproduce the measured (sub-) mm fluxes 
of 253-1536b. This is due to the fact that the emission of 
such small and dense disks becomes optically thick and, 
as a consequence, underestim ates the relativel y large 
(sub-)mm fluxes of 253-1536b (jTesti et al.ll200l[ ). Note 
that these possible values for the outer radius of the 253- 
1536b disk are all significantly lower than the estimated 
radiu s of the Roche lobe, i.e. >100 AU (from iPaczvnskil 
Il97ll using our estimates for the stellar masses, and the 
projected physical separation as a lower limit for the bi- 
nary semi-major axis). This means that the material in 
the disk lies well inside the stable zone in the Roche lobe. 

Another important parameter is the disk inclination, 
defined as the angle between the disk axis and the line-of- 
sight. By taking the ratio of the two projected disk axes 



in the HST images, we estimated an inclination i ~ 55° 
for the 253-1536a disk. This procedure cannot be ap- 
plied to the smaller 253-1536b disk, which has not been 
detected by HST. The X-Shooter spectra show that the 
253- 1536b star is significantly less extincted than its stel- 
lar companion. This indicates that the 253-1536b disk is 
likely less inclined than the companion. For this disk 
we considered a range of possible inclinations between 
0° and 50°. 

With the parameters outlined above we determined the 
spectral index j3 of the dust opacity coefficient by fitting 
the long-wave SEP with the two-layer disk models (see 
IRicci et all l2010af ) . Table [2] shows the constrained f3- 
values for the possible different R ou t and i for the 253- 
1536b disk, and in the p — 1 case. The difference between 
the estimated f3- values for the 253- 1536b disk is given by 
the different contribution of the inner optically thick re- 
gions to the total emission. The decrease of R ou t, or 
the increase of z, which increases the line-of-sight opti- 
cal depth of the disk, makes this contribution more im- 
portant. As a consequence, j3 becomes larger, i.e. the 
spectrum of the dust emissivity steepens, to compensate 
for the opposite effect given by optically thick emission. 
To better quantify the value of j3 for the 253- 1536b disk, 
high-angular resolution imaging is needed to directly con- 
strain its outer radius, and therefore the impact of the 
optically thick inner regions to its total emission. The 
important point to be noticed here is that for all the 
possible values of the disk parameters (this result is un- 
changed for possible other p- values between and 1.5), 
the j3- value constrained for 253- 1536b is larger than for 
253-1536a. 

The /3-index carries information on the size a max of 
the l argest grains in th e dust population of the outer disk 
(e.g. lNatta et al.l l2007). j3- values lower than about 1 — 1.5 
can only be explained with the presence of grains as large 
as at least 0.1 — 1 mm. Converting an estimate for f3 into 
one for a max is particularly difficult because of our igno- 
rance on the phsyical/chemical properties of the probed 
dust. However, for all the dust models considered in the 
literature, a general anticorrelation between j3 and a max 
is generally obtained for j3 < 1.5. Under the assumption 
that the chemistry /shape of the dust grains in the two 
disks is the same and considering that the dust evolution 
models in disks with the constrained physi cal properties 
predi ct a slope q which is nearly identical ([Birnstiel et al.l 
l20lll ). this means that the observational data indicate 
that the 253-1536a disk, with a lower /3, contains larger 
grains in its outer regions than 253 -1536b. By consider- 
ing the same dust model adopted in lRicci et al.l (|2010al lbL 
l201ll ). and varying the a ssumed g-value between 2 and 
3 (see IRicci et al.ll2010a[ ). the inferred a max in the 253- 
1536a disk is larger than in the companion disk by factors 
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TABLE 2 
Disk parameters in the p = 1 case. 



Source 


Rout 


i 


P 




,outer 


^"dust, outer 




(AU) 


(°) 




(g/cm 2 ) 


(K) 


253-1536a 


280 


55 


0.5±0.2 


1.5 


-4 


18 • 


-29 


253-1536b 


40 





1.0±0.2 


30 - 


-60 


18 


- 25 




40 


30 


1.4±0.2 


50 - 


- 100 


18 


- 25 




40 


50 














60 





0.7±0.2 


10 - 


-20 


15 


- 20 




60 


30 


0.8±0.2 


11 - 


- 22 


15 


- 20 




60 


50 


0.9±0.2 


12 - 


-24 


15 


-20 



Columns 5,6: ranges of gas surface densities and dust t emperatures constr ained in the disk outer regions (Section 14. 2 j) . These values have 
been obtained by adopting the dust model described in Ricci et al. ( 2010a). However, note that for our analysis, only the ratio of the 
surface densities and dust temperatures in the two disks is relevant. A model with R ou t = 40 AU and i = 50° cannot reproduce the 
observed fluxes for 253-1536b. 



of about 2 — 10, depending on the values for i? ou t an d i 
for the 253-1536b disk. 

4.2. Testing the models of dust evolution 

Recent physical models of dust evolution which include 
coagulation and fragmentation of dust grains predict a 
local relation 



(i) 



^max — 2 



between the maximum grain size allowed by fragmen- 
tation, the gas surface density E gas , the viscosity a v i SC - 
parameter ([Shakura fc Sunvaevlll973f ) , the thermal speed 
c s , the critical velocity Uf above which two grains frag- 
ment after colliding , and the grain density p s (see 
iBirnstiel etHI l2010al lbh . 

By assuming that the chemical composition and shape 
of dust grains in the two disks of the binary system is the 
same, the highly uncertain term uf / p s in Eq. [T] cancels 
out when taking the ratio of a max in the two disks: 



&max,l ^ ^gas,l C s,2 <^visc,2 ^ ^gas,l ^dust,2 <^visc,2 
&max,2 ^gas,2 C s C^ v isc,l ^gas,2 2~dust,l C^visc,! 



, (2) 



where T^ust is the dust temperature in the disk midplane. 
This means that observational constraints on the ratio 
of amax in different disks allow to test models of dust 
evolution in disks without basing the whole analysis on 
parameters whose values are very uncertain. 

Equations [T] and [2] are valid locally in the disk. How- 
ever, as noted above, our observations do not have 
enough angular resolution to properly resolve the disk 
structure and directly constrain the radial dependence 
of quantities like e.g. the dust surface density. The esti- 
mates obtained above on the dust properties in the two 
disks are referred to the disk regions which dominate the 
emission at long wavelengths, i.e. the outer regions. We 
can therefore attempt a comparison between the predic- 
tions of the dust evolution models and our observational 
results in the disk "outer regions" , which have to be prop- 
erly defined. 

Since the two disks have very different outer radii, the 
spatial regions probed by their (sub-) mm SED are differ- 
ent. According to the two-layer disk models, more than 
50% of the total sub-mm emission from 253-1536a comes 



from regions with stellocentric radii r > 100 AU, whereas 
for 253- 1536b the same fraction of emission comes from 
r > 20 - 30 AU. As shown in Table [2] for the p = 1 case, 
in these regions the 253-1536b disk is denser than 253- 
1536a by a factor of about 10 — 30, and colder by a factor 
of about 1.2-1.4. By including these values in Eq. [2j and 
assuming the same a v i SC - value in the two disks, it is ev- 
ident that the models of dust evolution would predict 
significantly larger grains in the 253- 1536b disk, which is 
contrary to what was derived in the last section. This 
result does not change even when considering other pos- 
sible p- values for the power-law index of the disk surface 
density between and 1.5. 

To reconcile the prediction of the models of dust evo- 
lution with the ratios of a max reported in Section 14. 1\ 
the viscosity a v i SC -parameter in the 253- 1536b disk has 
to be larger than in the companion disk by more than 
a factor of 10. Physically, this is because to explain 
the smaller grains observed in the outer regions of the 
253- 1536b disk, the turbulence velocity, which is roughly 
proportional to y/a v { sc , has to be very high to increase 
the relative velocities between solid particles thus making 
fragmentation more efficient. 

Although different a v i SC -values could in principle be 
present in the two disks, magneto-rotational simulations 
of protoplanetary disks predict that larger values of a v i sc 
are ty pically obtaine d in environments with lower den- 
sities ([Gammiel [l996) , which is the opposite of what re- 
quested by the dust evolution models to explain the ob- 
servational results for the 253-1536 binary system. The 
fact that larger grains are seen in environments with 
lower densities probably suggests that radial motion of 
particles in the disks, a phenomenon which is not in- 
cluded when deriving this prediction, plays a fundamen- 
tal role in the redistribution of s olid particles in proto - 
planetary disks. This is also what IBirnstiel et al.l (|2010bl ) 
suggested to reconcile the predictions of these same mod- 
els with the measured sub-mm fluxes of isolated disks in 
Taurus and Ophiuchus SFRs. 

High angular resolution and high sensitivity observa- 
tions with ALMA and EVLA will allow us to test the 
dust evolution models locally in the disk, and possibly 
to probe viscosity with high spectral and spatial resolu- 
tion observations of gas. 



Based on observations collected at the European 
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Fig. 1. — ONC map at 450 /im with SCUBA/JCMT ( Johnstone & Bally 1999). White crosses indicate the H ST- identified proplyds, 
whereas the star marks the location of 6 1 Ori C. The location of the 253-1536 binary system is shown. 
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Fig. 2. — Map of the 253-1536 binary system. In color the HST image in the Ra F658N filter ( R icci et al. 2008), whereas white contours 
represent the EVLA continuum data at 6.9 mm. Contour lines are drawn at 3, 6 and 9cr, where a = 0.035 mJy/beam. The white ellipse 
in the lower right corner indicates the size of the EVLA synthesized beam. 



